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The energy and centrality dependence of the charged multiphcity per participant nucleon is 
shown to be able to differentiate between final-state saturation and fixed scale pQCD models of 
initial entropy production in high-energy heavy-ion collisions. The energy dependence is shown to 
test the nuclear enhancement of the mini-jet component of the initial conditions, while the centrality 
dependence provides a key test of whether gluon saturation is reached at RHIC energies. HIJING 
model predicts that the rapidity density per participant increases with centrality, while the saturation 
model prediction is essentially independent of centrality. 
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Signals of the formation of a quark-gluon plasma in 
high-energy heavy-ion collisions are sensitive to the initial 
condition of the dense matter created in the early stage 
of the collisions. Global bulk observables at the Rela- 
tivistic Heavy-ion Collider (RHIC) such as the rapidity 
density of multiplicity and transverse energy provide im- 
portant constraints on those initial conditions. The first 
data on dNch/drj in Au-|-Au reactions at ^/s — 65 and 
130 AGeV were reported in Ref. We focus here on 
how the systematics (energy and centrality dependence) 
of this observable can be used to differentiate between 
competing entropy production mechanisms. The entropy 
is especially interesting because it remains approximately 
conserved during the complex dynamical evolution of the 
system if local equilibrium is maintained. In longitudi- 
nal boost invariant hydrodynamics the rapidity density 
of particles is in fact conserved. 

Theoretical estimates of the initial condition at RHIC 
vary over a wide range |^ because of the (as yet) 
unknown interplay between soft and hard mechanisms 
of multiparticle production in nuclear collisions. Phe- 
nomenological models ranging from soft string models 
§,|, perturbative QCD (pQCD) models |-|l§, and clas- 
sical Yang-Mills approaches have been proposed to 
predict the produced entropy. In these models there are 
basic physical parameters, such as the nuclear size de- 
pendence of the scale, po{^/s, A), separating perturbative 
QCD and nonperturbative components, that control the 
magnitude of entropy production. In multiple mini-jet 
pQCD-inspired models, such as HIJING, the unknown 
nuclear gluon shadowing and parton energy loss due to 
final state interactions also lead to a factor of ~ 2 vari- 
ation in the predictions for the final charged hadron ra- 
pidity density in the central region |]. With the first 
measurements of the dNch/drj from the PHOBOS exper- 
iment [Q and the other experiments at the Relativistic 
Heavy-ion Collider (RHIC), the theoretical uncertainties 
can soon be considerably reduced. 

The measured dNch/dri{\ri\ < 1) for central collisions 



was reported to agree within statistical and systematic 
errors with the default HIJING1.35 |^ predictions. Of 
course, one cannot conclude that the dilute HIJING1.35 
initial conditions (with dNgiue{pT > 2 GeV)/dy « 250 
at ^/s = 130 AGeV) are correct without much more ex- 
tensive differential experimental studies, especially of the 
shape of dNch / dr] and the high px hadron spectra ||l2|,|l3| . 
HIJING does not include final state rescattering except 
for schematic jet quenching. Other models, such as in 
EKRT jioj, that assume local equilibrium and hydrody- 
namic expansion can also reproduce the first data start- 
ing from much higher initial conditions. The aim of this 
letter is to emphasize that the centrality or impact pa- 
rameter dependence of the charged particle rapidity den- 
sity provides a key observable that, combined with other 
differential measurements, can significantly narrow the 
current band of uncertainty of the initial conditions pro- 
duced at RHIC and search for evidence of novel gluon 
saturation phenomena j^,|l^,|ll[ or dynamical screening 
effects |l^ at high density. 

In this study, we use HIJING Monte Carlo model ver- 
sion 1.35. This model incorporates pQCD to compute 
multiple mini-jet production and uses the Lund string 
model to describe soft beam jet fragmentation and 
hadronization of the jets. In default HIJING1.35, nu- 
clear shadowing of gluons is assumed to be identical to 
the observed quark shadowing. The large pt quarks are 
assumed to have energy loss of dE/dx = 1 GeV/fm, and 
jet quenching is modeled by a gluon splitting scheme. HI- 
JING unitarizes the minijet cross section via an eikonal 
approach p5[ 



NN 



= j dH[i 



(1) 



where ajet{s) is the inclusive jet production cross section 
with Pt > Po, (J soft accounts for soft interaction cross 
section and Tj^i\[{s,b) is the nucleon- nucleon geometri- 
cal overlap function. The two critical physical parame- 
ters of this model, po and a soft, are adjusted to fit the 
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measured cross sections and dNch/drj for inelastic pp{p) 
collisions at high energies. In versions 1.35 and below 
Duke-Owen (DO) parameterization of parton distri- 
butions is employed. If more recent parametrization of 
parton distributions in nucleons is used, one has to use 
an energy-dependent po in order to fit the cross section 
and hadron rapidity density in pp or pp. However, the 
results will remain the same as presented in this letter us- 
ing IIIJING1.35. The important point to keep in mind is 
that the cut-off scale po is fixed in HIJING by pp{p) data 
and is assumed to be independent of A or the centrality. 

Figure 1 compares dNch/drj per pair of participant nu- 
cleons in pp{pp) and central Au + Au collisions as func- 
tions of colliding energy. For nuclear collisions the two 
solid curves correspond to IIIJING1.35 calculations with 
(upper) and without (lower) jet quenching. We note that 
the effect of jet quenching on the total hadronic multi- 
plicity in HIJING model only becomes significant above 
> 100 GeV. 

From Fig. 1 both pp{p) and AA collisions appear 
to have a component that is approximately constant 
{n)soft ~ 1-3 plus a logarithmic energy dependent com- 
ponent. The constant component arises from the soft 
(low transverse momentum) hadron production due to 
beam jet string fragmentation in HIJING. This soft par- 
ticle production is also proportional to the number of par- 
ticipants, and therefore its contribution to dNch/drj per 
participant is independent of A in heavy- ion collisions. 
Including the contribution from minijet production on 
the other hand, the total hadronic rapidity density ac- 
quires the form 



dNa! 

drj 



= {Npart){n)soft + f {-^binary} 



is) 



(2) 



where aj'^{s) is the averaged inclusive jet cross section 
per NN collision in A A collisions and / w 1.2. We 
have checked that HIJING results without jet quench- 
ing indeed follows the above scaling. Since the de- 
fault HIJING includes nuclear shadowing of the gluon 
distribution, cr^'^ is generally smaller than cj^f . For 
central Au + Au{b = — 3fm) collisions the averaged 
number of participants is (Npart) ~ 350 and the aver- 
aged number of binary collisions per participant pair is 
'^.{Nbinary) / {Npart) ~ 4.7 at RHIC energies. Because of 
the s dependence on the minijet production cross section 
(T^f (s), the energy dependence of the hadronic multiplic- 
ity is amplified by '-^ A^/^ in heavy-ion collisions relative 
to pp{p) by the binary nature of semihard processes (see 
also Fig. 3 below). The rise of the multiplicity per partic- 
ipant relative to Super Proton Synchrotron (SPS) is con- 
sistent, within the large (mostly systematic) error bars, 
with the predicted binary scaling of the hard component 
from pp to AA via Eq. (0) . 
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Charged particle rapidity density per participating 
baryon pair versus the cm. energy. The PHOBOS data 
0] (filled circles) for the 6% most central Au-I-Au are 
compared to pp and pp data (open symbols) [p7|-p^ and 
the NA49 P6 -f P6 (central 5%) data (H] (filled square). 
HIJING1.35 (sohd) with (upper) and without jet quench- 
ing (lower) and EKRT (dot-dashed) predictions are also 
shown. 



Parton saturation can occur in both initial and final 
state interactions but they both give the same kind of A- 
dependence of particle production. The initial state sat- 
uration model |^,|llj is based on the nonlinear Yang-Mills 
field effect arising from gg g, in which case the satura- 
tion scale is determined by plat^^ = (9/16)C^7rQ!s AxG 
where Ca = 3, as is the strong coupling constant and xG 
is the gluon distribution per nucleon at x = 2psat/\/s 
and Q = Psat- Numerical solution of the Yang-Mills 
field illll could also include final state saturation and 

In this paper we will 



dynamical screening effects |14 



only compare our calculation with EKRT |10| model 
of final state saturation which gives a definite predic- 
tion of both energy and A dependence of rapidity den- 
sity of charged multiplicity. In this model, the pQCD 
growth of low pt gluons is cut off below a saturation 
scale, po{^fs,A) = p^at ~ 0.2AO i3(^s)O i9 at which 
dNg/dy=pl,,R\ 

Assuming direct proportionality between parton and 
the final hadron number, the energy and atomic number 
dependence of the total hadronic multiplicity per unit ra- 
pidity in central A + A collisions is estimated in EKRT 
ra as 



dNch 
dy 



(6 = 0) 



llAGA'^'^Ws) 



0.4 



(3) 



To compare to HIJING calculation and the PHOBOS 
data for pseudo-rapidity density, we scale the above re- 
sult by a factor 0.9. As shown in Fig. 1, the EKRT model 
estimate of the energy dependence of the entropy density 
is remarkably close to the HIJING results and is also 
consistent with the PHOBOS data. 
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However, a critical difference between these models is 
that the hadronic multiplicity per participant actually 
decreases in the EKRT model with the atomic number of 
the colliding nuclei due to the saturation requirement. It 
therefore considerably overestimates the hadronic multi- 
plicity if extrapolated down to pp{p). This is in contrast 
to the HIJING model of minijet production without sat- 
uration, where the multiplicity per participant increases 
with A according to Eq.(|^). 



Au+Au Collisions 




The charged particle (psuedo) rapidity density per par- 
ticipating baryon pair is shown as functions of the num- 
ber of participant baryons for central Au + Au(b < 3fm) 
colhsions at Vs = 56, 130 and 200 AGeV. Results of HI- 
JING1.35 (solid histograms) are compared to pp{p and 
PHOBOS data and to EKRT predictions (dot-dashed). 

Since extrapolating the EKRT model down to pp{p) 
collisions cannot be justified, it is more useful to study 
centrality dependence of heavy- ion collisions. We cal- 
culate next the centrality dependence of the hadronic 
multiplicity per participant in order to emphasize the 
difference between fixed scale and saturation models of 
entropy production. Shown in Fig. 2 are the dNch/drj per 
participant pair as functions of (Npart) at three different 
energies together with the RHIC data by PHOBOS ex- 
periment iQl and the pp{p) data. The HIJING results in- 
crease monotonically with the number of binary collisions 
per participant {Ntinary) / (Npart) in an intuitive way as 
given in Eq.(|^). The slope increases with energy since 
the hadronic multiplicity is proportional to jet cross sec- 
tion which has a significant energy dependence. For very 
peripheral collisions, the results agree with pp{p) data be- 
cause that is how the model parameters of HIJING are 
constrained. Naively, {Nunary) / {Npart) ^ (Npart)^^^- 
The deviation from such a simple dependence in HIJING 
calculation is due to a combined effect of jet quench- 
ing and dilute edges in Wood-Saxon nuclear distributions 
used in HIJING. 

The dot-dashed lines are the predictions extrapolated 



from EKRT saturation model Eq. (||), assuming the de- 
pendence on atomic number in central collisions roughly 
the same as the number of participant pairs Npart I for 
a fixed A. A generalization of the EKRT approach to 
model local saturation is expected not to change our esti- 
mate of its centrality dependence qualitatively. Contrary 
to HIJING predictions, the saturation model gives in- 
creasing multiplicity per participant toward more periph- 
eral collisions. While the extrapolation to the highest im- 
pact parameter (low participant number) domain is du- 
bious, a general feature of saturation models is expected 
to be a weakly decreasing or constant dependence on cen- 
trality in semi-peripheral to central collisions. The up- 
coming experimental data should easily distinguish these 
two widely different predictions. 

We emphasize the power of the combined energy and 
centrality dependence in Fig. 2 to search for novel nu- 
clear saturation effects. This also applies at the Large 
Hadron Collider (LHC) energies. The higher the energy 
the greater the difference between the predicted central- 
ity dependencies. A threshold for saturation can be iden- 
tified experimentally by looking for a region with nega- 
tive derivatives d{dNch / drjNpart) / dNpart < above some 
Npart > Ncrit{\/s). Ordinary dynamical Screening cffccts 
on the other hand, could be recognized by a gradual re- 
duction of the positive derivatives predicted by HIJING 
with increasing centrality. 

Of course more differential data will eventually provide 
much stronger constraints on models. Measurements of 
large px hadron spectra can provide for example infor- 
mation related to jet quenching |p^ , |l3| . Another impor- 
tant differential observable is the rapidity dependence of 
dNch/drj/ (Npart) ■ The rapidity dependence reflects the 
x-dependence of the gluon distribution function in nu- 
clei and should be strikingly different in saturation mod- 
els. In Fig. 3, the rapidity dependence of this quantity 
is shown for the HIJING model. Since most of the gluon 
production is from small x region, one can see that the 
nuclear enhancement in the hadronic rapidity density per 
participant is mainly restricted to the central rapidity re- 
gion in this model. The width of this region is roughly 
determined by the scale pq as Ay ~ Iii{^/s/pq). If satura- 
tion occurs, a different rapidity dependence of the nuclear 
enhancement may result. The nuclear size and energy de- 
pendence of the shape of that distribution as compared 
to different model calculations adds an important observ- 
able in the search for possible saturation effects 

Finally we note that the study of the centrality depen- 
dence of nuclear enhancement per participant depends 
on the experimental ability to measure or deduce Npart- 
At present this is done assuming that dNch/drj ex Npart 
in the fragmentation regions 77 > 3. In Fig. 3, we see that 
at least in the HIJING model the nonlinear nuclear de- 
pendence is indeed confined to 77 < 3. However, detailed 
studies of the centrality dependence of correlations be- 
tween dNch/drj at = and 77 > could help reduce 
systematic errors. 
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The charged particle rapidity density per participating 
baryon pair as functions of rapidity as predicted by HI- 
JING for central Au + Au{b < 3fm) (solid) and pp 
(dashed) coUisions at = 56, 130 and 200 AGeV. 
Also shown are experimental data for pp collisions at 
= 200 GeV. 



In summary, we have shown how the energy and cen- 
trality dependence of the hadronic multiplicity densities 
in the central region of high-energy heavy-ion collisions 
can be used to constrain the mechanisms responsible for 
producing the initial conditions in such reactions. The 
energy dependence of the charged hadron multiplicity has 
been shown by PHOBOS to be enhanced relative to pp 
reactions and is consistent with the onset of pQCD dy- 
namics driven by the binary nature of semi-hard (mini- 
jet) processes. Though similar energy dependence is pre- 
dicted by very different pQCD based models, the central- 
ity dependence can differentiate between them. In prin- 
ciple, parton saturation is expected to occur [p| jlo|Jlll ] at 
asymptotic high-energy in collisions of very heavy nuclei. 
The interesting question is where this occurs in practice. 
We proposed that the centrality or A dependence of the 
hadron multiplicity at RHIC and higher energies may 
help to answer this question. 

Note added in proof: Since the submission of this let- 
ter, the first data |^l| on the centrality dependence of the 
multiplicity has become available. Other calculations, 
motivated by this work, have also appeared, e.g., p2|,p3|. 
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